Study on Unconventional Superconductors via Angle-resolved Specific Heat 
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The gap function in unconventional superconductors may vanish at points or lines in momentum 
space, permitting electronic excitations, termed nodal quasiparticles, to exist at temperatures well 
below the superconducting transition. In the vortex phase, the presence of nodal quasiparticles 
should be directly observable through the variation of the heat capacity with the angle between a 
magnetic field and the location of the zeroes of the gap. The heat capacity of candidate non-magnetic 
unconventional superconductors Lu(Y)Ni2B2C were found to exhibit fourfold oscillations with field 
angle, the first such observation. The observed angular variations are in quantitative agreement 
with theory, confirming that quasiparticles are created via Doppler shifts at nodes along < 100 >. 
Anomalous disorder effects have been also observed in the field-angle dependent heat capacity C p (a). 
In a slightly disordered sample, anomalous secondary minima along < 110 > appeared for fioH > 
1 T, leading to an eightfold pattern. The coexistence of an anisotropic superconducting gap and 
nonlocal effects is shown to drive the anomalous behavior. These results demonstrate that field- 
angle-dependent heat capacity can be a powerful tool in probing the momentum-space gap structure 
in unconventional superconductors such as high T c cuprates, heavy fermions, borocarbides, etc. 



I. INTRODUCTION 

Most superconductors behave conventionally; elec- 
tronic excitations are suppressed by the BCS gap causing 
the electronic heat capacity, for example, to be exponen- 
tially small at temperatures well below the superconduct- 
ing transition. Recently, a number of superconductors, 
cupratesi and heavy-fermion metals^ among them, have 
been found to be unconventional in that they exhibit gap- 
zero (or nodal) points or lines in momentum space. Elec- 
tronic excitations - nodal quasiparticles (nqp) - are then 
observed at low temperatures, giving rise to power-law 
rather than exponential, behavior. Unlike gapless super- 
conductivity, which can occur in conventional supercon- 
ductors, the Fermi momenta of these quasiparticles are 
restricted to nodal regions of the Fermi surface, giving a 
strong directional dependence to various physical prop- 
erties. 

In 1992, Yip and Sauls proposed a nonlinear Meiss- 
ner effect in the penetration depth of cuprate high-T c 
superconductors^ Since the effect uses the field-direction 
dependence of the supercurrent in locating the positions 
of the nodal lines (or points) of an unconventional gap in 
momentum space, it is a stronger test of the symmetry 
of the superconducting state than power-law behavior for 
the penetration depth. The search for the fundamental 
manifestation of d-wave symmetry has been very contro- 
versial from an experimental point of viewm£&Z& Li et 
al. questioned the observability of the nonlinear Meissner 
effect. 9 They introduced nonlocal electrodynamics for a 
d-wave superconductor and found that the local nonlin- 
ear Meissner effect is not observable for fields below a 
crossover scale H* describing the competition between 
nonlinear and nonlocal effects in the Meissner state. For 



most orientations of the screening current, the crossover 
is of the same order as or greater than the lower critical 
field H cl . 

In conventional type II superconductors, the vortex 
cores play a crucial role in transport properties because 
a relatively high density of bound states is created there 
by the suppression of the local order parameter and ex- 
tended quasiparticle states are completely gapped. In 
a new symmetry class such as cuprate superconductors, 
however, Hirschfeld et al. argued that extended quasi- 
particle states are dominant over vortex bound states 
because of the presence of the nodes^ The anisotropy 
in the gap in momentum space fixed to the crystal 
axes will induce an angular anisotropy in the current re- 
sponse due to the coupling between quasiparticles and 
supercurrent, leading to the fourfold angular variation 
in field-angle thermal conductivity of unconventional su- 
perconductors with d-wave order parameter. In 1995, Yu 
and Salamon et al. demonstrated for the first time that 
the heat transport of YBCO oscillates with in-plane field 
angleiii*i£ Subsequent experiments have confirmed the 
feature and thermal conductivity has been established as 
a tool that can probe the momentum-specific gap nature 
in unconventional superconductorsiiiiiiiii Such measure- 
ments, however, are complicated by the competition be- 
tween Andreev-scattering and Doppler-shift effects. 

A much more direct, and far more difficult, exper- 
iment is to detect the low-energy density of states of 
the nodal quasiparticles through its modulation by an 
in-plane field. Heat capacity provides such a measure. 
Working in the 2D limit, Vekhter et al. have shown that 
the density of states of a d-wave superconductor exhibits 
four-fold oscillation with field angle with respect to crys- 
tal axes. 16 From the experimental stand point, the search 
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for the oscillation has been puzzling piLiSi The expected 
30 % of oscillation amplitude is within the reach of exper- 
imental error, but the effect has proven elusive. Whelan 
and Carbotte argued that the Zeeman energy splitting 
fiH competes against the anisotropyiiS At a critical field 
H c , Zeeman and Doppler energy scales become compa- 
rable and the DOS variation with field angle vanishes. 
Recently Won and Maki used a semi-classical approxi- 
mation and extended the 2D Vekhter model to a 3D su- 
perconductor with modulated cylindrical Fermi surface 
by including quasiparticles with an out-of-plane momen- 
tum component The 3D effect causes a strong suppres- 
sion of the oscillation amplitude to around 6 %, making 
the observation much more difficult. 

In this paper, we review the direct observation of a phe- 
nomenon that should be common to many so-called un- 
conventional superconductors, those whose Cooper pairs 
have symmetries more complex than originally envisaged 
by Bardeen, Cooper and Shrieffer in their classic the- 
ory of superconductivity. If, in particular, the pairs have 
the symmetry similar to d-orbitals in atoms, there will be 
certain directions in momentum space in which electronic 
excitations, termed quasiparticles, can arise without an 
energy gap. 

In section 2, the experimental method, sample prepa- 
ration and characterization will be briefly explained. In 
section 3, the field- angle heat capacity of LUN12B2C is 
shown to be similar to results reported previously for 
YNI2B2C1SI The magnetic field dependence will be an- 
alyzed by a Won-Maki extension of the Vekhter et al. 
calculation. The heat capacity dependence on magnetic 
field angle will be discussed in terms of 3D nodal quasi- 
particle theory and used as a piece of confirming evidence 
that LUM2B2C and YNi2B2C belong to a class of super- 
conductors with nodes in gap function. An eightfold pat- 
tern appeared in a slightly disordered Lul22I supports 
the coexistence of nonlocal effects and anisotropic gap 
effects. Then, concluding remarks will follow in section 
4. 



II. EXPERIMENTS 

The electronic specific heat is a fundamental physical 
quantity that measures the electronic density of states 
(DOS) directly and is a bulk property. As noted in the 
Introduction, it was recently suggested that the specific 
heat can be used to study the angle-resolved gap struc- 
ture of unconventional superconductors by studying the 
DOS variation with magnetic field angled Considerable 
efforts have been expended to detect the field-angle vari- 
ation arising from the anisotropic gap structure in un- 
conventional superconductors, especially in YBCOii2*i& 
Those efforts, however, have suffered from intrinsic ex- 
perimental limitations. First, a very limited set of field 
directions were investigated. Typically, the crystal-axes 
and diagonal directions were chosen because they give 
the largest contrast in DOS for pure d x i_ y i or d xy pairing 
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FIG. 1: Sample assembly in ac calorimetry 



symmetry. If the pairing symmetry is the mixture of two 
order parameters, such as d+s— wave, then the two direc- 
tions do not represent gap-maximum or gap-zero (nodes) 
directions because the nodes move away from the diago- 
nal directions. Second, temperature-variable heat capac- 
ity intrinsically contains large background contributions 
such as lattice vibrations or thermally excited quasipar- 
ticles. The large backgrounds cause an uncertainty in 
extracting the field- induced part, and can easily obscure 
the small angular variation. 

With those problems in mind, we built a new specific 
heat probe designed exclusively to study the gap struc- 
ture of superconductors in momentum space. Instead 
of varying temperature as was done in many previous 
works, we fixed temperature and magnetic field intensity 
and made the field angle as an only variable. This ex- 
perimental scheme enabled us to study the small angular 
oscillation, which is typically less than I % of the mea- 
sured heat capacity. Some of the key ingredients to this 
scheme are as following. First, temperature and mag- 
netic field should be very stable. Otherwise, even a small 
fluctuation could obscure the small change with angle. 
A Physical Properties Measurement System (PPMS) by 
Quantum Design was used as a cryogenics platform to 
satisfy those conditions. Second, a stepping motor with 
a gear ratio 1:141 was used to rotate the sample and was 
controlled by a driver from PPMS. Instead of measuring 
limited directions, we were able to study the angle de- 
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FIG. 2: Upper panel shows the field-angle heat capacity of 
optimally doped YBCO at 1.7 K and 7 T. The phase shift 
with field angle is shown in the lower panel. In both panels, 
the squares describe the data measured by the first pair of 
thermocouple wires and the circles by the second pair. 



pendence of heat capacity in increments of 3° or smaller. 
Third, the heat capacity measurement should be sensitive 
enough to pick up the minute change with field angle. In 
addition, the measurement must be able to handle small 
single crystals because their weight is usually less than 
1 mg. AC calorimetry was chosen because it has a sen- 
sitivity of up to 0.01 % of the measured heat capacity 
as well as a capability to measure even microgram (/ig) 
samples i^i 

The ac method used here is an adaptation of the Sul- 
livan and Seidel technique^ Fig. 1 shows the sample as- 
sembly The front face of the sample was coated with col- 
loidal graphite suspension (DAG) thinned with isopropyl 
alcohol to prevent a possible change of the optical absorp- 
tions of the sample. The sample was weakly coupled to 
the heat bath through helium gas and suspending ther- 
mocouple wires. At high temperatures, the thermal cou- 
pling can be adjusted mainly by controlling the amount of 
the helium gas surrounding the sample, while the thermal 
leak through thermocouple wires should be taken into ac- 
count at cryogenic temperatures. As a heating source, in- 
frared laser light (A = 789.2 nm) was used. The light was 
modulated electronically to make a square-wave pulse 
and was guided into the PPMS sample space through 
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FIG. 3: The averaged field-angle heat capacity of YBCO at 
1.7 K and 7 T. The solid line shows a theoretical calculation 
of fourfold oscillation with the oscillation amplitude of 3.4 % 
(see text). 



an optical fiber. The oscillating heat input incurred a 
steady temperature offset (or dc offset) against the heat 
bath with an oscillating temperature superposed. The 
measured oscillation temperature T ac was converted to 
heat capacity by the relationship C oc 1 /T ac in a proper 
frequency range. 

In our field-angle heat capacity measurements, two 
thermocouple pairs were used to measure the temper- 
ature oscillations T ac simultaneously. The two sets of 
Au/Fe and chromel thermocouple wires were configured 
so that they were 90° apart, which enabled us to rule out 
the field-angle contribution from the thermocouple wires 
due to field-dependent thcrmopower of Au/Fe. If the an- 
gular oscillations come from our experimental setup, the 
two signals will be out of phase while they will be in 
phase with each other if they come from sample. As an 
example, we studied an optimally doped YBCO (Fig. 2). 
The squares in the top panel describe the angular heat 
capacity measured by the first pair of the thermocouple 
wires and the circles by the second pair. The two sets of 
data were measured simultaneously at 1.7 K in 7 T with 
magnetic field rotating within the Cu02 plane of YBCO. 
Both of them show clear two-fold oscillations with field 
angle measured against a-axis. Note that the oscillations 
are out of phase with each other, indicating that they 
come from the thermocouple wires. The bottom panel 
of Fig. 2 shows the phase shift with field angle. There is 
no correlation with the heat capacity shown in the upper 
panel, confirming that the two-fold oscillation of the heat 
capacity is not related to intrinsic sample properties. 23 

In Fig. 3, we show the averaged heat capacity of YBCO 
(2C ave = C\ + C2). As expected, the two-fold compo- 
nents are almost suppressed. Now that the artifact from 
the experimental setup is removed, we can study the in- 
trinsic sample property. The solid line describes a fit by 
a cusped fourfold variation with oscillation amplitude of 



3.4 %, i.e., C(a) = C + c 4 (l + T| sin2a|). Here, a is the 
magnetic field angle with respect to the a-axis and T is 
the oscillation amplitude. The periodicity is 7r/4 and the 
minima are located at (110) directions, which are consis- 
tent with d x 2_ y 2— wave order parameter in YBCO where 
the nodes are located along (110). The field-angle oscil- 
lation, however, is obscured by noise which is comparable 
to the angular change. 



III. OBSERVATION OF NODAL 
QUASIPARTICLES AND NONLOCAL EFFECTS 
IN THE NONMAGNETIC SUPERCONDUCTORS 
LU(Y)NI 2 B 2 C 

The discovery of RNi 2 B 2 C (R = Y, Lu, Tm, Er, 
Ho, and Dy) captured attention of the superconductiv- 
ity community because of its relatively high transition 
temperature (T c ) and the coexistence of magnetism and 
superconductivityi2ii25i2& The highest T c 's in the boro- 
carbide family are 16.5 K and 15.5 K in the nonmagnetic 
members with Lu and Y elements respectively. Amid 
the controversy over the pairing symmetry of the boro- 
carbides, there is growing evidence that the gap function 
is highly anisotropic and possibly has nodes where the 
gap becomes zero in momentum space. In specific heat 
measurements, a power law behavior was observed in the 
temperature dependence and the electronic coefficient 
'y(H), to the extent it can be extracted, follows a square- 
root field dependencei2L2& Thermal conductivity mea- 
surements showed low energy excitations at as low as 70 
mK and an anomalous magnetic field dependence similar 
to the unconventional superconductor UPt3 but differ- 
ent from the exponential dependence characteristic of an 
isotropic s-wave gap^S Compelling evidence for the pres- 
ence of nodes along (100) directions has been reported 
from both z-axis field-angle thermal conductivity^ and 
field- angle heat capacity measurements of YNi 2 B 2 CiSi 

Recently, a magnetic field-driven flux line lattice (FLL) 
transition has been observed both in the magnetic mem- 
ber (Er) and in the nonmagnetic members (Y and Lu) 
of the tetragonal borocarbidesi2ii22i2iiM The transition 
from square to hexagonal vortex lattice occurs due to the 
competition between sources of anisotropy and vortex- 
vortex interactions. The repulsive nature of the vor- 
tex interaction favors the hexagonal Abrikosov lattice, 
whose vortex spacing is larger than that of a square lat- 
tice. The competing anisotropy, which favors a square 
lattice, can be due to lattice effects (fourfold Fermi sur- 
face anisotropy)^ or unconventional superconducting or- 
der parameter^ In this section, we present evidence that 
both the gap anisotropy and nonlocal effects coexist and 
the coexistence is the origin of the anomalous crossover 
from a fourfold pattern to an eightfold pattern in C(a) 
with increasing field in LuNi 2 B 2 C. 
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FIG. 4: Temperature dependence of the in-plane resistivity 
of sample A (crosses) and sample N (circles) . 



A. Disorder effects 

The best samples of LuNi 2 B 2 C exhibit behavior iden- 
tical to that reported earlier for YNi 2 B 2 C?22i but disor- 
dered samples are significantly different. Three single 
crystals of LuNi 2 B 2 C from different batches were stud- 
ied. Sample N was annealed at T = 1000 C° for 100 
h under high vacuum, while sample A and C were not 
annealed. The annealed crystal, sample N, showed the 
highest T c (16.1 K) and residual resistivity ratio (RRR) 
of 27, while the unannealed sample A recorded the low- 
est T c of 15.5 K and RRR of 16 (see Fig. 4), suggesting a 
judicious postgrowth annealing may improve the sample 
quality. The resistivity at T c is 2.34 and 1.44 /Ltfi-cm for 
samples A and N, corresponding to mean-free paths of 
144.5 and 234 A respectively. Assuming that 16.1 K is 
the T c for a pure sample, sample A is equivalent to 0.8 % 
of Co doping on the Ni site; i.e., Lu(Nii_ a: Co a: ) 2 B 2 C with 
x = 0.008.21 Sample C with T c = 15.8 K is less disordered 
than sample A. 

For an isotropic s-wave superconductor, the contribu- 
tion to the density of states (DOS) comes from the lo- 
calized states in the vortex core. Since the density of 
vortices is proportional to the magnetic field, the field- 
induced DOS is linearly dependent on the field and (3=1. 
For superconductors with line nodes, delocalized states 
leaking through the nodes dominate the DOS. When 
we integrate the contribution in a vortex unit cell and 
multiply by the density of vortices, the electronic DOS 
is proportional to in an intermediate region, i.e. 
H c i <C H <C H C 2, where H c2 is upper critical field. This 
so called Volovik effect gives (3 = 0.5 for line node and 
j3 = 1 for point node>2& Recently, there have been several 
reports that s-wave superconductors such as CeRu2>22> 
and NbSe2^£ behave like d-wave superconductors, show- 
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ing sub-linear magnetic field dependence in specific heat. 
Ichioka et alM> extended the Volovik calculation to the 
entire region, i.e. H c i < H < H C 2- While the Volovik ef- 
fect considered the delocalized states only, they included 
the localized states bound in a vortex core and the quasi- 
particle transfers between vortices. In their numerical 
calculation, the DOS dependence in the d-wave pairing 
states gives an exponent (3 = 0.41, rather than Volovik's 
P = 0.5. As for those with anisotropic s-wave pairing 
symmetry, the exponent becomes 0.67, which is clearly 
distinguishable from 0.5 or 0.41. 

The disorder effects are also manifested in the mag- 
netic field dependence of the electronic specific heat C e ; . 
In pure samples (sample C and N), the low-temperature 
specific heat shows a \f~H dependence in the mixed state, 
but there is a deviation from the \ftl dependence at a fi- 
nite field in the disordered sample A. Fig. 5(a) shows the 
heat capacity of sample C at 2.5 K with increasing (cir- 
cles) and decreasing (crosses) magnetic field along [100] 
direction. The dashed line which represents the least 
square fit of Co + b(H — Hq)@ best describes the field 
dependence with (3 — 0.46. Here Co is the zero-field heat 
capacity and the offset Hq is essentially a lower critical 
field (H c i) and 0.1 T was used for the least-square fit. 

Figures 5(b)-(d) show the heat capacity of the disor- 
dered sample A at 2.5, 4, and 8 K respectively. Low-field 
data were best described by the square-root field depen- 
dence (dashed lines) and the arrows indicate the points 
where the data deviate from the fit. The deviation field 
shows a systematic increase with temperature, i.e, 0.8 T 
at 2.5 K, 1.8 T at 4 K, and no clear deviation at 8 K. 
Above the deviation field, the data fall below the VH 
line. 
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FIG. 5: (a) Heat capacity of sample C at 2.5 K for H ||[110]. 
(b)-(d) Heat capacity of sample A at 2.5, 4, and 8 K, respec- 
tively for H || [110]. Dashed lines are H 1 ^ 2 fits and arrows 
indicate the deviation points from the H 1 ^ 2 dependence. 38 



B. Angle-resolved Specific Heat 



In addition to the v H dependence in unconventional 
superconductors, the density of states N(H,a) also de- 
pends on the orientation of the superfluid flow with re- 
spect to the nodes, which results in angular oscillation in 
C(a)iiSili The supercurrent flow around a vortex leads 
to a Doppler energy shift, SE = v s • Kk-F , where v s is the 
velocity of the superfluid and Kkp is the Fermi momen- 
tum of nodal quasiparticles. When the field direction is 
normal to the plane containing nodes, the DOS is the av- 
erage over the whole Fermi surface, leading to a square 
root field dependence. When the field is in the nodal 
plane, however, the Doppler shift has a field-direction 
dependence as well, SE w ^ sin {3 sin(</> — a)w& Here <f> 
is an azimuthal angle of the gap node and (3 is a vortex 
current winding angle. The energy scale associated with 
the Doppler effect is defined as 



Et, 



ahv* 



$ 



1/2 



(1) 



where the geometrical constant a is order of unity, v* = 
yjv a v c is the average Fermi velocity, and $0 is the flux 



quantum. The variable p — r/R, where r is the distance 
from the vortex core and 2R is the intervortex distance. 
Working in the 2D limit, Vekhter et al. calculated the 
DOS , N fa (Ni + N 2 )/2, for an in-plane magnetic field 
when uj,Eh <C Ao where to is the Matsubara frequency, 
and for four nodes at angles a n from orthogonal axes in 
a plane : 



Nj(u), h, a) 
~N~o 



— 1 

Ej 



1 

2X 1 



[(1 + 2x 2 ) arcsinx + 3a; vl 



{x = uj/Ei > 1) 

5 (x<l), 
(2) 



where i=l,2, E\ — Eh\ sin(a„ — a)\ and E2 — 
Eh,\cos(a n — a)\. The energy E\ is for nodes close to 
the field direction, while E2 is for nodes far from it. A*o 
is the normal-state density of states. 

When the magnetic field is along a nodal direction, 
a = a n + 7nii / 2 where m is an integer; the two nodal po- 
sitions perpendicular to the field contribute fully, while 
the two parallel to the field do not contribute at all be- 
cause the Doppler shift is an inner product between the 
superfluid velocity v s and the quasiparticle momentum 
Kkp. When the field is along a maximal gap direction 
(a = a n + (2m + 1)tt/4), in contrast, all four nodal points 
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contribute equally to the DOS with a factor of y2 less 
than the full contribution. When we add all the contribu- 
tions, the total DOS has minima along nodal directions 
and maxima along maximal gap directions. At T = K, 
the oscillation is cusped with a contrast of 1/ Vz between 
the minima and the maxima and is independent of the 
magnetic field. At a finite temperature, however, the 
sharp contrast is washed out due to thermal effects, lead- 
ing to more rounded minima and the DOS oscillation now 
depends on the magnetic field. 

The oscillation amplitude depends on the gap geom- 
etry and dimensionality of the superconductors. The 
Vekhter et al. model above assumes the Fermi momen- 
tum of nodal quasiparticles (nqp's) to be restricted to 
the nodal plane. A more realistic calculation, therefore, 
is to allow nqp Fermi momentum out of the nodal planes, 
which decreases the amplitude of the angular variation. 20 
We account for the 3D effect by replacing the Doppler- 
related energies Ei and E 2 by Ef D = Eh [sin 2 (a n — a) + 
cos 2 9} 1 / 2 and E\ D = E h [cos 2 (a n — a) + cos 2 9} 1 / 2 respec- 
tively and integrate the DOS over polar angle 9: 

1 f 2 " 1 

N(w,h,a) = — / -[N 1 (w,h,a,9)+N 2 (w,h,a,6)]d6. 

(3) 

The qualitative features are the same as those in the 2D 
case except that the oscillation amplitude decreases from 
41 % to 6 % in 3D superconductors at T = K. The 
3D effect was indeed confirmed experimentally in the 
field- angle heat capacity measurement of YNi2B2C»22i We 
note that a similar effect is predicted for an (s+g)-wave 
superconductor^ 

The upper panel of Fig. 6 shows the field-directional 
angular dependence of the total heat capacity at 2 K 
and 0.6 T, where the transverse magnetic field was ro- 
tated within the basal plane of LUM2B2C. The total 
heat capacity consists of constant, 2-fold, and 4-fold con- 
tributions: C to tai{a,H) = C + C 2 {a,H) + C 4 {a,H). 
The field-independent constant Co is due to nonmag- 
netic contributions, such as the lattice heat capacity and 
thermally excited nqp's, and is determined experimen- 
tally from C(T) at T. The 2-fold contribution C 2 (a, H) 
comes from our experimental setup and has a functional 
form of c 2 cos 2a. The Au/Fe and chromel thermocouple 
wire is a major source of this contribution, but misalign- 
ment of the basal plane of the sample against the field di- 
rection could also lead to the 2-fold component because of 
the anisotropy between in-field ab-plane and c-axis heat 
capacities. The dashed line in the upper panel of the 
Fig. 6 is C 2 (a, H); the 2-fold signal is about 40 % of the 
4-fold component at 0.6 T and increases with magnetic 
field. 

The lower panel of the Fig. 6 shows the 0.6 T data af- 
ter subtraction of the 2-fold and the lattice components 
: C 4 (a,H) = C tota i(a,H) - C 2 (a,H) - C . The 4-fold 
variation is clearly seen. To make a quantitative analysis, 
we fit our data to C 4 (a,i?) = 04(H) (1 +T\ sin2a|). The 
coefficient c^(H) and the angular contrast T were treated 
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FIG. 6: Field directional dependence of the heat capacity of 
sample A at 2 K and 0.6 T. The field angle a is measured 
with respect to the a-axis. The top panel shows the total 
heat capacity (open circles) and the 2-fold component, C2(a) 
relative to baseline of 4.9 (dashed line). The bottom panel 
shows the same data after subtraction of the 2-fold and lattice 
contributions, C4 = C to tai — C2(a) — Co- The data were 
best described by a cusped function, Ci(a,H) — C4(H)(1 + 
T\ sin2a|) with c 4 = 1.45 and T = 0.0415. 



as field dependent fitting parameters. The solid line in 
the lower panel of Fig. 6 shows the fit with T — 0.0415 
at 0.6 T. The sharp minima along < 100 > indicate 
that there exist gap minima or nodal structures along 
those directions. The nodal positions and the small os- 
cillation amplitudes are consistent with previous reports 
of YNi2B2di^°. The field-angle heat capacities of sam- 
ple C and N at low fields are essentially the same as that 
of sample A. 

Figure 7 gives 3D surface plots of the field-angle heat 
capacity of sample C at 2.5 K. The x-axis is the in-plane 
field direction a with respect to the crystal axes and the 
y-axis is the applied magnetic field intensity in tesla. The 
z-axis is the four-fold heat capacity normalized by the 
Doppler-induced heat capacity: Ci(H,a)lci(H). Fig- 
ure 7(a) is a numerical calculation of the 3D nodal quasi- 
particlc theory 16,20 and Fig. 7(b) are the experimental re- 
sults of sample C. In the quasiparticlc theory, the only ad- 
justable parameter is the Doppler-related energy Eh(v*) 
(see Eq. 1). In the numerical calculation, we used a = 1, 
the gap maximum A = 1.76fcsT c , and the average Fermi 
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FIG. 7: The 3D surface plots describe magnetic field-angle (x- 
axis) and field-intensity (y-axis) dependence of the in-plane 
specific heat ratio at 2.5 K: Ci(H, a) /04(H). Fig. (a) is a 
numerical calculation of the specific heat due to nodal quasi- 
particles in 3D system and Fig. (b) is the experimental data 
of the single crystal sample C. The experimental data were 
FFT low-pass filtered with a cut-off frequency of 0.03 Hz with 
angles being considered as time. Both plots (a) and (b) were 
constructed with 0.5, 1, 1.5, and 2 T data sets. 



velocity v* — y/v a v c = 1.3 x 10 7 cm/sec. The extracted 
Fermi velocity is similar to that from the band struc- 
ture calculation (^Jv a v c ~ 1.2 x 10 7 cm/sec)^ The 3D 
plots of sample N show similar behavior to that of sam- 
ple C requiring a larger Fermi velocity, i.e., v* = 1.5 x 10 7 
cm/sec^ The larger value may be due to the uncertainty 
in determining the field-induced heat capacity. 

In field-angle heat capacity measurements, there are 
three relevant energy scales, which are the Doppler en- 
ergy scale Eh, the paramagnetic energy fj,H, and the ther- 
mal energy of ksT. In the mixed state, the paramagnetic 



1.55 




angle 

FIG. 8: Angle-resolved heat capacity of sample A at 2 K and 
a 4 + 4 pattern (solid line). The magnetic fields are 0.6, 1, 
1.5, and 2 T, from top to bottom panel. 



energy is of order 1 K and the Doppler energy scale is 
~ 10 K. The latter is much higher than the experimental 
temperature 2 K. Therefore, the Doppler effect on the 
nodal quasiparticles is dominant over the thermal and 
the paramagnetic effects in this regime, indicating that 
the 3D nodal quasiparticle theory of the Doppler effect is 
applicable as was evident from the reasonable agreement 
between experiments and the numerical calculations. 

Figure 8 shows the heat capacity of sample A as a 
function of field angle at 2 K and 0.6, 1.0, 1.5, and 2.0 T 
from top to bottom panel. At 1 T, surprisingly, addi- 
tional minima develop along < 110 >, producing two 
sets of fourfold patterns or 8-fold, an effect not observed 
in either sample C or N with higher T c 's. The crossover 
field from the fourfold to the (4 + 4) pattern of sample A 
lies between 0.6 and 1 T, which is also the point where 
the heat capacity of sample A deviates from the y/~H de- 
pendence (see Fig. 5). With increasing field, the splitting 
gradually disappears and the field-angle heat capacity re- 
covers its fourfold pattern above 4 T. We also measured 
the field-angle heat capacity of sample A at 4 K to check 
if the anomalous peak splitting persists at higher tem- 
peratures (Fig. 9). The fourfold pattern now persists to 
1 T, evolving into two sets of fourfold patterns above 2 T. 
The 4 T data at 4 K have a shape similar to the 2 T data 
at 2 K. The crossover field H s i increases with increasing 
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FIG. 9: Angle-resolved heat capacity of sample A at 4 K and 
a 4 + 4 pattern (solid line) . The applied magnetic fields are 
1, 2, 3, and 4 T, from top to bottom panel. 



temperature. 

It is worth noting that the anomalous 8-fold pattern 
occurs only at sample A, which has half the electronic 
mean free path of the sample N, while the T c is slightly 
decreased. According to the nonlocal theory by Kogan 
et oi.; !5 the hexagonal- to- square FLL transition depends 
on the electronic mean free path I and the superconduct- 
ing coherence length £. Gammel et al. found that a mere 
9% of Co doping onto the Ni site in Lul221 can make 
the FLL transition field 20 times higher than that of the 
pure matrix for H | [001]. 4 » The FLL transition field 
for pure Lul221 is around 0.2 Tj 4 ^ which is well below 
our measurement range, suggesting the nonlocal effects 
would not influence the field-angle heat capacity of purer 
samples. In contrast, the disorder in sample A increases 
the transition field to at least twice higher than that of 
sample i.e., to a field relevant in the field-angle heat 
capacity measurement of sample A. 

When the magnetic field is rotated within the a&-plane, 
the transition field may differ with different field direc- 
tions because of the different nonlocal range, i.e., 
The two different transition fields can be characterized 
by H s i for < 110 > and H S 2 for < 100 >. As a magnetic 
field rotates within the a&-planc for H s \ < H < H S 2 
, the FLL will experience a structural change (or dis- 
tortion), i.e. hexagonal for H | [100] and square for 
H || [110]. Since the borocarbides have nodes on the 
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FIG. 10: Oscillation amplitude F of sample C (circles) and 
sample N (squares). The dashed line is from the 3D nodal 
quasiparticle theory with vf = 1.5 x 10 7 cm/sec. The bottom 
panel contrasts the oscillation amplitude due to the nonlocal 
effects (circles) and gap anisotropy (squares). 38 



Fermi surface, the DOS will differ depending on the FLL 
structured The negative deviation from the vH above 
H s i in the heat capacity of sample A (see Fig. 5b) also in- 
dicates that the DOS of the hexagonal FLL is larger than 
that of the square FLL. In addition to the anisotropic 
gap modulation^ the DOS will oscillate due to the FLL 
anisotropy, leading to a (4+4)-fold pattern. 

When H < H s \ or H > H S 2, the field- induced pa- 
rameter C4 can be considered as constant and the DOS 
oscillation is just from the gap anisotropy. When H s \ < 
H < H S 2, in contrast, C4 is not independent of field di- 
rection any more, but oscillates due to the anisotropic 
FLL transition field: C4 = pl(l + j\sin2(a — 45)|). C4 
has maxima along < 100 > and minima along < 110 >. 
Then, the field-angle heat capacity can be written as 

C 4 (a) = pl(l + y\sin2(a - 45)|)(1 + r|sm2a|), (4) 

where pi is a field-dependent fitting parameter. The 
value r represents the oscillation due to gap anisotropy 
in pure samples (see Fig. 10). The nonlocal effects give 
rise to a 45°-shifted 4-fold pattern and are accounted for 
by 7. The solid lines in Fig. 8 and 9 are least square 
fits of Eq. (4) and represent the data very well. The 
oscillations due to the nonlocal effects (7) and the gap 



9 



8 



' 1 1 1 

LJ 


' I ' 

Sample A 


n r.. ■ ■ 




Normal State | 


£^ <110> - T 
<100> • 

h r" 

s1 


■ 

■ 

■ 

■ 

■ 



4 8 12 16 

T(K) 



FIG. 11: H — T phase diagram of sample A with lowest T c . 
Hsi is the FLL transition for H ||[110] and H s2 for H || [100] . 
The dotted lines are guides to the eyes. Triangles and squares 
were sketched to show the corresponding FLL. 38 

anisotropy (T) at 2 K are compared as a function of mag- 
netic field at the bottom panel of Fig. 10. The FLL effect 
7 increases sharply above 0.6 T and decreases gradually 
to zero at 4 T, indicating that the low field corresponds 
to H s i and the high field to H S 2- It is interesting to 
note that Miranovic et al. came to a similar conclusion 
by solving quasi-classical Eilcnberger equation for nodal 
superconductors with anisotropic Fermi velocity^ 

The H — T phase diagram of the disordered sample A 
is shown in Fig. 11. H s i is the FLL transition field for 
H || [110] and H s2 , for H || [100]. The increase in the 
transition fields with increasing temperature is consistent 
with the nonlocal prediction where the nonlocal range de- 
creases with raising Ti2& The difference in the FLL tran- 
sition fields for H || [100] and H || [110] can be due to 
the upper critical field H C 2 anisotropy between the two 
directions. Since the effective nonlocal range is propor- 
tional to the coherence length £, the H c i anisotropy leads 
to a higher FLL transition field for H | [100] than that 
for H || [110]. The FLL transition occurs when the in- 
tervortex distance is comparable to the effective nonlocal 
radius. The observed ratio H s2 /H s i w 4 is larger than 
the predicted ratio 2^ but is smaller than the reported 
ratio of 10 in YN^L^C. 3 ^ The difference between exper- 
iments and theory may attest to the need that we take 
into account both the nonlocality and the anisotropic gap 
nature of the borocarbides. 



IV. CONCLUDING REMARKS 

We reviewed the first direct evidence for a variation 
of the DOS of nodal quasiparticles (nqp's) of uncon- 
ventional superconductors Lu(Y)Ni2B2C. The four-fold 



field-angle oscillation in the heat capacity of the boro- 
carbides is due to the nqp's experiencing a field-induced 
Doppler energy shift. The 3D superconductivity of the 
borocarbides is reflected in the small oscillation ampli- 
tude of 4 %, not 40 % expected for 2D superconductors. 
A dramatic change has also been observed in a slightly 
disordered Lul221. A fourfold pattern in C(a) changed 
to an eightfold or two sets of fourfold pattern at the FLL 
transition field, evidencing that both nonlocal effects and 
an anisotropic superconducting gap coexist in the boro- 
carbides. 

Even though it is clear that there are nodes along 
< 100 > directions in the borocarbides, the type of 
nodes has yet to be answered. Based on a theoretical 
calculation, 4 ^ it is asserted that the rapid decrease of 
the oscillation amplitude of the angle-dependent ther- 
mal conductivity with polar angle 8 is due to a point 
node gap located along [100] ~ Miranovic et al, how- 
ever, claimed that the oscillation amplitude decreases 
sharply with 9 for both vertical line node and the s + g 
point node cases, suggesting that the conical field rota- 
tion experiments may not be able to distinguish the two 
cases. 45 In order to determine the precise gap structure, 
instead, they suggested the field dependence of the os- 
cillation amplitudes at low fields. The oscillation mono- 
tonically decreases with increasing field for point nodes, 
while it has a maximum for line nodes. We plan to mea- 
sure conical field-angle heat capacity and extend our mea- 
surements to lower fields, which could shed more light 
on the controversy of the nature of nodes. A system- 
atic study on disorder effects through Co doping on Ni 
site, LuN^-zCozL^C, will clearly elucidate the correla- 
tion between nonlocal effects and anisotropic gap effects. 

It is gratifying that other field-angle heat capacity have 
been followed after our pathbrcaking worki^S Deguchi 
et al. applied the technique to the spin-triplet super- 
conductor Sr2Ru04 and found a fourfold oscillation for 
H _L c— axis, providing decisive information on the multi- 
band superconductivity and gap structures^ Aoki et al. 
have studied the heavy-fermion superconductor CeCoIn5 
and found a clear fourfold oscillation in C(a) with min- 
ima along < 100 >, suggesting the presence of line nodes 
at those directions^ Their interpretation of the super- 
conducting gap symmetry as d xy — wave, however, con- 
tradicts that (d x 2_ y 2) by thermal conductivity data^S 
It is interesting to note that the angle-resolved specific 
heat directly measures the zero-energy density of states 
(ZEDOS). On the other hand, the thermal conductivity 
data k(o) necessarily involve both the ZEDOS and the 
quasiparticle scattering time r(a)fii making it difficult 
to identify the gap node direction in some cases. 

Finally, we mention that the angle-resolved specific 
heat measurement opened a new venue to better under- 
standing unconventional superconductivity. Even though 
only a few classes of unconventional superconductors 
have been studied by this technique j2£L224L22i new fea- 
tures on their superconducting properties have been re- 
vealed, indicating a whole new aspect of superconductiv- 
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ity may await to be found. 
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